On 
On 

00 



in 



c 
o 
o 



> 

oo 



O 



T3 

c 
o 
o 



X 



Toward a Theory of Orbiton Dispersion in LaMnO,3 

Vasili Perebeinos and Philip B. Allen 
Department of Physics and Astronomy, State University of New York, Stony Brook, New York 11794-3800 

(February 1, 2008) 

At 750K, LaMnC>3 has a cooperative Jahn- Teller (JT) distortion, with Mn atoms in distorted 
oxygen octahedra. This lifts the degeneracy of the singly-occupied e g orbitals of the Mn 3+ ions, 
which then become orbitally ordered. We use a minimal model to describe the ordered phase at 
T — 0. The on-site Coulomb repulsion U is set to infinity. There are two electronic orbitals and three 
oxygen vibrational coordinates per unit cell. In addition to spin excitations and phonons, the model 
has electronic excitations consisting of mis-orienting orbitals on Mn ions. Neglecting coupling to the 
oxygen displacements, the gap to such excitations is 2A = lQg 2 /K where g is the electron-phonon 
coupling and K is the oxygen spring constant. When static oxygen displacements are coupled, this 
excitation becomes a self-trapped exciton with energy A, half the JT gap. Adding dynamic oxygen 
displacements in one-phonon approximation introduces dispersion to both the (previously Einstein- 
like) phonons and the orbital excitons ("orbitons"). One of the phonon branches has zero frequency 
at k = (0,0,7r). This is the Goldstone mode of the JT broken symmetry. 

I. INTRODUCTION 

The fascinating properties of Manganite perovskites have recently been rediscovered [0-0- LaMnOs has Mn 3+ 
ions with electronic configuration 3d 4 in the high-spin t\ q +e\+ state. The single electron in the doublet e g level is 
Jahn- Teller (JT) active. A cubic to orthorhombic distortion B occurs at Tjt=750K which lifts the degeneracy with 
alternating sign in alternating cells, and drives orbital ordering 0]: x and y-oriented e g orbitals alternate in the 
x — y plane. This cooperative JT transition yields a layer structure which then determines at the much lower Neel 
temperature (Tn=240K) an antiferromagnetic "A-type" (AFA) |j| spin-ordering with ferromagnetic spin order in the 
x — y planes. We believe that the principal term in the Hamiltonian which drives the cooperative JT transition is 
coupling of Mn orbitals to motions of oxygen atoms along the directions of the (approximately 180°) Mn-O-Mn bonds. 
In a previous paper M we have discussed a model for the JT transition and used this model both to describe the 
T = phase of pure LaMnOs and the formation of small polarons when hole-doped as in Lai-ajCa^MnOa. Our model 
is related to the Rice-Sneddon S model for BaBiOs, and was introduced for LaMnOs by Millis M. 




FIG. 1. Schematic of Jahn- Teller distorted LaMnOs x — y plane with orbitals ordered 



Our calculations are done at T=0 and x=0. Previous authors flld-h2] have studied orbital dynamics using models 
which omit JT and electron-phonon effects, focussing on dynamics induced by hopping. This leads to complicated 
orbital dynamics strongly dependent on spin order because of Hubbard and Hund energies. Our model emphasizes a 
very different kind of physics. Because we take U = oo and x—0, there is no hopping and the spins are decoupled from 



other degrees of freedom. The remaining electron-phonon physics is still very rich, and we think it will be dominant 
at small x. Other treatments of the small x orbital dynamics exist (see ref. |13| and papers cited there), but the large 
ratio Tjt/2"n argues against omitting the electron-phonon physics at small x. 




FIG. 2. Tetragonal Brillouin zone of the Jahn- Teller distorted phase of LaMnOa 



II. MODEL AND GROUND STATE 



Each cubic cell has 5 degrees of freedom, two electronic and 3 lattice dynamical. The electronic degrees of freedom are 
occupancies of the Mn e g orbitals. In Van Vleck notation [Q ip2 and ^3 are the familiar x 2 — y 2 and 3z 2 — r 2 functions. 
These can be rotated by 45° to an alternative orthogonal basis, the symmetrical pair ipx, ^Py = (—^3 ± V^Vv^- 
These orbitals point in x and y directions respectively, with a moderately strong lobe in the z direction. 

Lattice dynamical degrees of freedom are labeled by u{£ ± x), u(£ ± y), and u{£ ± z), which refer to oxygen 
displacements along the x, y, and z directions of the oxygen atoms half-way between the Mn atom at site £ and the 
next atom in the ±x, ±y, and ±z directions, respectively. The notation u(£ — x) is equivalent to u(£' + x) for the Mn 
site £' , which is translated from £ by one unit in the — x-direction. 

It is convenient to use both the Cartesian labeling given above, and a Van Vleck labeling: 



^Qi^) = [u{£ + x) - u{£ - x)] + [u{£ + y)~ u(£ - y)] + [u{£ + z) - u{£ - z)} 
Q 2 (£) = [u{£ + x) - u{£ - x)] - [u{£ + y)- u(£ - y)] 
V3Q 3 (£) = 2[u{£ + z) - u{£ - z)] - [u{£ + x) - u(£ - x)\ - [u{£ + y) - u{£ - y)} 



(1) 



Our starting Hamiltonian TL has a lattice term 7Yl and an electron-phonon term TL op . The lattice term does not have 
a useful expression in the Van Vleck notation, and is best written as 



Wl 



E 



f^ + A) + J^ + A) 



(2) 



On the other hand, the electron-phonon term is best written in Van Vleck operators, 



Wep = -<?E( C xW,4W 



Q 2 (£) -Q 3 (£) 
-Q s (£) -Q 2 (£) 



cx(£) 
cy{£) 



V2gJ2Qi( £ ) ^ x {£)c x {£) + c\.{£)c Y {£)) 



In the undoped case, each Mn site is singly occupied; we stay in the subspace where the operator relation 

c x (£)c x (£)+c f Y (£)c Y (£) = l 
holds. The degenerate e g orbitals are defined to have zero energy. 



(3) 



(4) 



In adiabatic approximation, the ground state of this Hamiltonian has a doubled unit cell and a net average oxygen 
displacement (Q 2 (£),Q3(P)) of the type Q(cos8,sm8) exp(i<f ■ I). The wave-vector q — (tt,tt,tt) is favored, but the 
angle 8 in the (Q 2 , Q3)-plane is undetermined. The observed distortion, of the Q2-type (8 = 0), is fixed by terms 
left out of our Hamiltonian. Once a Q2~type order parameter is established, the wave-vector q can equally well be 
(it, 7r, 0), and this is what is observed. The lattice is now bi-partite. We label sites by A when exp(itf • £) is 1, and by 
B when it is -1. According to Eqn. pi the A sites will be occupied by X orbitals and the B sites by Y orbitals. The 
alternating occupancy of these orbitals on A and B sublattices is shown schematically in Fig. ([j]) . 

Oxygen atoms move off-center by ±2g/K, giving a JT distortion Q 2 = 8g/Kexp(iq- £). This lowers the energy 
of each occupied orbital by 8g 2 /K = A and raises the empty orbitals by the same amount. The distortion raises 
the lattice energy by 4Ng 2 /K, which is half the energy lowering of 8Ng 2 /K from J aim- Teller energy. The lowest 
electronic excited states ("orbitons") are misoriented orbitals which cost 16g 2 /K = 2A «2eV. 

III. ORBITON HAMILTONIAN 

We are interested in the lowest excited states of the model. We allow vibrations of oxygen atoms around their 
distorted positions. Therefore we are no longer making the adiabatic approximation. Lattice vibrations will renor- 
malize both the ground state and the orbiton excitations. By flipping an orbital we create an orbital defect. We 
can define a creation operator a^(£) for the orbital defect on the A site as c Y (£)cx(£), and similarly on the B site 
tf(£) = c x {£)cy{£). One can show that these operators obey Fermi commutation relations in the subspace of singly 
occupied sites. Furthermore one can show that on the A sublatticc number the operator of the misoriented orbitals 
c y (£)cy(£) is the same as a)(£)a{£); similarly on the B sublattice c\A£)cx(£) is the same as b^(£)b(£). 

Using these definitions we can rewrite the Hamiltonian Eqns. 0, in terms of the lattice fluctuations u(£ + A) 
defined as uq(£ + A) + Su(£ + A), where uq(£ + A) = ±2g/K in x or y direction and in the z direction. 



4Ng 2 /K + H cp +H h = 2 A ( ^ a^(£)a(£) + ^ tf (£)b(£) j 



v 2 v ' 2M £ 

£,X—x,y,z 

2g (y,5Q 2 {£)a\£)a{£) - ^ 5Q 2 (£)b\£)b{£) j 

■ 9 ( E S Q»W ( at W + ffl W) + E WW ( &t W + 6 W) J ( 5 ) 

Wa ten ) 



The next step is into Fourier space. Define Fourier variables by 



°w - \/lE 



N 

k 



a k e ik -* 



fo w = /|E^' 



k-l 
Uk<i 

k 



M^A) = ^l£^^e- (6) 

fe 

where 5uk x = fk x + f—k 5 if ^ is on the A sublattice, and otherwise 5uk x = 5fc A + g_^ ■ The Fourier-transformed 
Hamiltonian has the form 

H = H° orh + H ph + H 3 (5Q 3 ) + H 2 (6Q 2 ). (7) 

The first two terms are unperturbed dispersionless orbitons and phonons with energy 2A and hio. 



H° rb = 2A J2( a t a k + b kb> 



H ph = hu Yl (ftj^+ald^) ( 8 ) 



The last two terms give the coupling of orbiton density to quantum fluctuations of the lattice of the <5Q3-type and 
<$Q2-type. 

n 3 (5Q 3 ) = K3 J2 [(4 + o-fc) (Mk) + 4("*0) + (b{ + b- k ) (A 2 {k) + Al(-k)) 

k 

n 2 (SQ 2 ) = k 2 J2 [( B ^ k ) + B I(-*)) p-( k ) + ( B ^ k ) + b|(-a)) p fc (fc)] (9) 

k 

We are now explicitly using the (ttttO) wavevector for Jahn- Teller oxygen distortions. New operators A, B, are 
introduced to represent the Fourier transformed SQ 3 (k) and SQ 2 (k) variables, and new operators p a and p b are used 
for orbiton densities on A and B sites: 

A 1 (k) = 2f kz (1 - e- lk >) - f kx - f ky + g k „e- ik * + g ky e~^ 



A 2 (k) -. 


= 2g kz (l - e~ ik ") - g kx - g ky + f kus e 


~ lkx +fk y e lk 


Bi(k) -- 


= fk x ~ fk y ~ 9k x e- ik " + g ky e- lk * 




B 2 (k) -- 


= 9k x - 9k y - .fk x e~ lk * + f ky e- tk y 




Pa(k) -- 


= VivE a U« 9 

q 




Pb{k) -- 







(10) 

1 /9 

The term due to SQ 3 couples orbitons to the phonon subspace with coupling constant k 3 — g (H/6Mcj) ' = 

1 /2 

(7iu>A/48) . Since the phonon energy is much smaller than the orbiton energy (hcu <C A) this term can be treated 
perturbatively. The term due to Q 2 couples phonon vibrations to orbiton densities with the coupling constant 
K 2 = 2g(h/2Muj) 1 / 2 — (TicjA/4) 1 / 2 . This term describes oxygen vibrations around new positions which relax by an 
amount proportional to the orbiton density p in order to lower the energy of the orbiton excitation. To deal with 
term we make a displaced oscillator transformation. In place of the phonon operators f kx , g kx defined above, we will 
use the following operators: 

F kx =fk x + ^(pa(-k)+ Pb (-k)e ik *) 

F ky =f ky -^(p a (-k) + p b (-ky k y) 

Gk * = 9k *~it M- fc ) eifc * + ^-*)) 

Gk y = g ky + £ (Pa(-k)e ik » + Pb (-k)) 

Fk z = }kz 

Gk> = g kz (ii) 

These new operators are bosonic. They were chosen to complete squares in the 7ip h +7i2 (SQ 2 ) part of the Hamiltonian, 
which becomes 



H° ph + H 2 (6Q 2 ) = hu ]T (FlF kk +G\G 

k,X=x,y,z 



^ E [(Pa(k) + Pb (k)e-^) x (KC) - (p a (k)e-^ + p h (k)) x (KC)] (12) 



where HC stands for Hermitean conjugate. 

We now truncate the Hilbert space to the subspace with only one orbiton present. Then the product of two adjacent 
orbiton annihilation operators gives zero. Then Eqn. U2 can be simplified further by noting that ^2 k [p a (k)p a (— k) + 
Pb{k)pb(—k)] is the same as J2 k [ at< (k)a(k) + tf (k)b(k)]; cross terms are zero in our subspace. The second term of Eqn. 
|l2| then combines with Tt® lb Eqn. (0) to renormalize the orbiton energy by — 4(n 2 ) 2 / hui = —A. The orbiton energy is 
now exactly half of the bare energy 2 A. The reason is obvious in a site representation. Flipping an orbital creates a 



local region where an orbital is surrounded in the xy plane by other orbitals of the same orientation. Oxygen atoms 
adjacent to the central (flipped) orbital will now prefer to locate in undistorted positions. This means that rather 
than raising the local energy from —A to +A with the flip, we have instead raised it from —A to 0. This is the sole 
influence of <5Q2-type fluctuations in the one-orbiton subspace. 

Now we transform the rest of the Hamiltonian to the displaced oscillator operators. Substituting new phonon 
operators ( |ll|) into H 3 (SQ 3 ) generates the same term as in (g), but with old phonon operators replaced with new 
ones, plus additional terms coupling orbiton operators to orbiton density operators. The new Hamiltonian has the 
following form: 



H = Tiorb + Hph + 7^1+7^2+7^3 

where the first two terms are the new bare orbitons and phonons 

Horb = A ^2 ( a l a k + b\b k \ 

k 

H ph = hw ]T (Fl x F kx +GlG kx 

k,X—x,y,z 

There are now three coupling terms left, 

Hi = Ks Z H Ai ( fc ) + b i A ^ k ) + hc ) 

k 

H 2 = n 3 J2 (<*Ui(- fc ) + & 14(- fc ) + HC 

k 

H 3 = n 2 Z ( cos kx ~ cos k y^ \ a kPb(-k) + b\p a (-h) + HC 



(13) 



(14) 



(15) 



Here the A and B operators are the same as defined in (M) after substition of the new phonon operators F and G in 
place of / and g. The new coupling constant k^ is ^K-iKz/Tw = A/VT2. The perturbation part of the Hamiltonian 
Hi + 7^2 + 7^3 couples orbitons and phonons to states with energies lying above or below by the order of A. These 
terms will be handled in the next section. 



IV. LOW ENERGY EFFECTIVE HAMILTONIAN 



In the Hamiltonian ([hi 15) orbitons and phonons are coupled directly through the Hi term. The term H3 does 
not couple to the phonons, but couples the one orbiton subspace to the two orbiton subspace. The remaining TL2 
piece couples the one orbiton and one phonon states to states with one more orbiton and one more phonon. Thus the 
Hamiltonian matrix acting on the one phonon and one orbiton subspace has the structure: 



H^ 



(hull Hi H 2a 
Hi AI 


\ 
H3 H 2 b 




( * (lph) 
*(lorb) 


H 2a (A + 2Tiuj) I 
H 3 
\ H 2b 



2AI 

(2A + Tiuj) I J 




*(lorb- 

*(2orb) 

V *(2orb- 



\ 

2ph) 
lph) J 



(16) 



The operator H 2 has matrix elements 7^2a when it couples the 1 phonon to the 1 orbiton plus 2 phonon subspace, and 
matrix elements H 2 b when it couples the 1 orbiton to the 2 orbiton plus 1 phonon subspace. The operator Hi couples 
1 orbiton to 1 phonon. It also couples 1 orbiton plus 2 phonons to 2 orbitons plus 1 phonon. However, this block does 
not influence the 1 phonon or 1 orbiton answers except in higher order, so these elements of the Hamiltonian matrix 
are omitted. 

This Hamiltonian can be perturbatively reduced to two separate low-energy effective Hamiltonians, one for single 
orbitons and one for single phonons: 



HI 



H 



Hi{E-huj) 1 Hi+H 2h {E-{2A + hcj)) 1 H 2h +H 3 {E-2A) l H 3 



Kh = ^Ph +Hi(E- Ay 1 Hi + ft 2a (E-(A + 2hu)) l H 2 



(17) 



The 2x2 effective Hamiltonian 7i°f b acts on wavefunctions built from two types of orbitons (al , 6j) | 0,0 > . Similarly 
the 6x6 effective Hamiltonian 7ip b acts on six-component wavefunctions (Fj A ,Gj A ) | 0,0 >, where A = x,y,z. The 
unperturbed ground state | 0,0 > with zero orbitons and phonons couples through the Ti.2 term to states of one 
phonon and one orbiton. This leads to a ground energy shift SEqs —< 0, | H.2 (E — (A + hui))~ H.2 | 0, >. 

The problem is now reduced to computation and diagonalization of a 2 x 2 and a 6 x 6 matrix. Note that we could 
have taken into account the 7ii part of the perturbation Hamiltonian by leaving it as an off-diagonal coupling term, 
giving an 8 x 8 effective Hamiltonian to compute and diagonalize. The error we make by folding it down as in l\li\ } 
has the same order of magnitude (fioj / A) 2 as the other errors we make in our perturbation theory. 

Our next task is to evaluate the matrix elements of the effective Hamiltonians (|l^) . After straight-forward algebra 
the result (with the ground energy shift SEgs = — 12k§AT/ (A + two) subtracted out) is: 

K% = (a - { -?$p\ i+ (^- + aT^)( (12 " 8 cos kz) f ~ 2 (cos kx + cos ky) &x 

/eff _ t, ,t , ( K 3 , K 3 \ ( Wo, 7i b 



Where a x is a Pauli matrix. The 3x3 matrixes 7i a and Tib are 

2 i + e *(fe. E -fe y ) -2 + 2e~ ife * 

H a = I l + e ^ k "- k ^ 2 -2 + 2e~ ik * 

-2 + 2e ife - -2 + 2e lfc - 8-8cos/c x 

-2cosfc 2: _ e -*fcy _ e «x 2e lk * - 2e l{k ^ k ^ 

H b =\ - e ^y- e - ik - -2cosk y 2e ik * - 2e l(k «- k > ] ) (19) 

2 e ~ ifc « _ 2e i ( kz ~ k *) 2e~ iky — 2e i ( kz ~ ky '> 

By solving ( jig ) for eigenvalues and eigenvectors we find the dispersion law of orbitons and phonons. In our 
approximation it is consistent to neglect Tluj relative to A in the denominators of (|l|). Then we obtain two dispersive 
modes for orbitons and two for phonons, while 4 phonon modes remains degenerate with eigenvalues E A '°' = hu. 
The dispersive modes are: 



-^orb " 


ii A hu , 

= -2 A+ U^ 


— 4 cos k z + cos k x + cos k y ) 


-C'orb " 


11 fwj , 
--2 A+ U^ 


— 4 cos k z — cos k x — cos k y ) 


-^ph = 


- Tiu) (6 — 

12 V 


4 cos k z + cos k x + cos k y ) 



E ph = Twj — — (6 — 4 cos k z — cos k x — cos k y ) (20) 

The Jahn- Teller distortion Q2 with a (tt, tt, 0) wavevector creates a simple tetragonal Brillouin zone, show in Fig. 
y. Fig. y shows the orbiton and phonon eigenvalues along principle symmetry lines. An interesting result is that 
the energy of one of the phonon branches goes to zero at wavevector k = (0, 0, tt). The ground state energy of our 
Hamiltonian is degenerate in (Q2, Q3) plane. It depends only on magnitude of the vector Q — (Q2, Q3), but not on its 
direction. We have chosen a Q2-type distortion as is observed in experiment. In our model, an infinitesimal distortion 
of the Q3 type does not cost any energy. This Goldstone mode occurs necessarily at the Z point of the Brillouin zone 
because only a q — (tt, tt, tt) order parameter is independent of the direction of ordering in (Q2, Q3)-space. 
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11/12 
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A T R W X 



FIG. 3. Dispersion relation of orbitons (upper two curves) and phonons of Jahn- Teller distorted LaMnOa. Energies are 
measured in units of the JT gap A. Four phonon branches decouple from orbitons and retain their bare energies hui, independent 
of wavevector. The parameter r\ = hui/A has been set to 0.25 for display purposes. A likely value of this parameter is 0.06. 

V. CONCLUSION 

In this work we discuss the dispersion of the orbiton excitation caused by coupling to phonons. We solve a simple 
but reasonably realistic model, which has previously been used MM to describe the insulating nature and magnetic 
ordering of lightly doped LaMn03. The solution makes a single occupancy approximation, allowing only one electron 
per Mn site. We restrict consideration to creation of only a single orbiton defect. We show that the electron-phonon 
interaction reduces the energy of the excitation by a factor of two, which is due to lattice relaxation. We also obtain 
a zero energy excitation or Goldstone mode which is due to a broken continuous symmetry of the Hamiltonian. We 
ignored two other sources of dispersion of the orbiton - the hopping terms (t 2 /U), and direct quadrupolar Coulomb 
coupling. These effects would probably cause a similar size dispersion to the one we considered, but would not have 
the major effect of reducing the JT gap by a factor of two which occurs because of coupling to phonons. 
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